s- Wave-Like excitation in the superconducting state of electron-doped cuprates with 

c!-wave pairing 
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An intrinsic physical mechanism, based on the doping evolution of the Fermi surface (FS), is 
explored to reconcile the contradictory experimental results on the superconducting (SC) pairing 
symmetry in electron-doped cuprates. It is argued that the FS pocket around (7r/2, 7r/2) has not yet 
formed until doping reaches about the optimal value. Therefore, in the underdoped regime, even if 
the SC order parameter is d wave which vanishes along the line k x = k v , the quasiparticle excitation 
gap is still finite and looks s-wave-like due to the absence of the FS across that line. This makes 
it possible, with d-wave SC pairing, to understand those experiments which evidenced the s-wave 
quasiparticle excitation. An explicit theory with consideration of both antiferromagnetic and SC 
orders is implemented to exhibit the FS evolution from underdoping to overdoping and associated 
with it the variations of the quasiparticle property, electronic density of states, and low temperature 
dependences of the physical quantities heat capacity and superfluid density. 

PACS numbers: 74.72. Jt, 74.20.-z, 74.25.Nf, 74.25.Bt 



I. INTRODUCTION 



The electron-doped (e-doped) cuprate supercon- 
ductors, typically Nd 2 - K Ce x Cu0 4 (NCCO) and 
Pr2_ K Ce :E Cu04 (PCCO), are an important component 
in the study of high-T c superconductivity. Compared 
to hole-doped cuprates, the e-doped ones have a much 
more robust antiferromagnetic (AF) phase and a much 
narrower superconducting (SC) phase. Moreover, there 
is strong evidence for the coexistence of both AF and SC 
orders around optimal doping x = 0.15^*^ although this 
issue is still under debated The transport measurements 
on PCCO^ suggest a quantum phase transition from 
antifcrromagnetism to paramagnetism at x ~ 0.17. The 
schematic phase diagram is plotted in Fig.^a). 

As the fundamental topic, the SC pairing symmetry 
has been intensely studied. While the s-wave symme- 
try was supported by the measurements like penetra- 
tion depthSil and tunneling^*2ii& the d-wave symmetry 
was suggested by the seemingly more extensive measure- 
ments including phase-sensitive Josephson junctions^ii^ 
angle-resolved photoemission spectroscopy (ARPES), 13 
specific heatji^ and also penetration depth>iSii& In addi- 
tion, there is other evidence to show the so-called non- 
monotonic ci-wave beh.aviosi2ii& or a transition from ol- 
io s-wave symmetryiSiSS with increasing doping. Thus 
the experimental data are rather controversial to support 
the s-wave or essentially (i-wave (no matter monotonic or 
not) pairing. In particular, a puzzle arised in the pen- 
etration depth measurements that different results were 
obtained by different groups and even the same group. 
For example, a finite SC gap was inferred by Kim et 
al£ over a wide range of doping levels from underdop- 
ing to overdoping, even though a transition from d- to 
s-wave behavior with increasing doping was claimed ear- 
lier by the same groupi^S On the other hand, completely 
contrary result, i.e., d-wave behavior at various doping 



levels, was reported by Snezhko et al^£ Why are the ex- 
perimental data so contradictory? Perhaps it is partly 
due to differences between fabricated samples. However, 
in this paper we intend to explore an intrinsic physical 
mechanism, based on the doping evolution of the Fermi 
surface (FS), to understand the diverse experimental re- 
sults. With it the apparent discrepancies might be largely 
reconciled. 
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FIG. 1: (a) The schematic phase diagram for electron-doped 
cuprates Nd^r^-xCexCuO^ The AF order persists beyond 



optimal doping x = 0.15. 



(b) The calculated order pa- 



rameters m ("x") and A ("o") at T = for three doping 
levels x = 0.13, 0.15 and 0.17 from the phenomenological t- 
t'-t"-U-V model to reproduce the experimental features. The 
lines are only for guiding the eyes. 



To expound the mechanism, we first examine the FS 
structure of e-doped cuprates. The ARPES measure- 
ments on NCC02i*2£ have revealed the intriguing FS 
evolution with doping. It was found that at low doping 
x = 0.04 only a small FS pocket appears around (it, 0). 
With increasing doping new spectral intensity around 
(ir/2, 7r/2) begins to emerge and becomes enough strong 
at optimal doping x = 0.15, indicating the formation of 
the new FS pocket around (7r/2,7r/2). The ARPES data 
have been explained within the two-band model in con- 
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sideration of the AF order^ 2 ^^ or by numerically con- 
sidering the AF short-range correlation. 26 According to 
the two-band analysis, at low doping only the upper AF 
band is crossed by the Fermi level around (tt, 0), leading 
to the small FS pocket around (w, 0). With increasing 
doping the AF order weakens so that the lower band, 
with the maximum at (7r/2,7r/2), approaches the Fermi 
level and contributes the new spectral intensity around 
(-7r/2,7r/2). The question is at which doping the lower 
band will be crossed by the Fermi level. As seen from the 
ARPES data on NCCO at x = 0.13 (Fig. 2 of Ref.|H, 
the lower band is still away from the Fermi level, although 
it is close to the latter and thus contributes weak spectral 
intensity around (ir/2,n/2) in the intensity map (Fig. 1 
of Ref. 22). In view that the lower band lifts continuously 
with increasing doping, we naturally imagine that it does 
not touch the Fermi level until x ~ 0.15. That is to say, 
in nearly the whole underdoped regime (x < 0.15), the 
real FS pocket around (tt/2, tt/2) has not yet formed. 

On the other hand, as far as we know, many experi- 
mental measurements to study the SC pairing symmetry 
only give direct information about the quasiparticle exci- 
tation gap, and not the SC order parameter. For exam- 
ple, from the low temperature dependence of the penetra- 
tion depth (exponential function or power law), one can 
only infer whether the quasiparticle excitation has a finite 
gap or not. It must be noted that a finite quasiparticle 
excitation gap, which was actually obtained by e.g. Kim 
et aZ.ji does not necessarily mean a nodeless SC order 
parameter. Here our essential idea is: The quasiparticle 
excitation gap may be finite and look s-wave-like, but 
the SC order parameter is d-wave and has nodes. This 
possibility becomes true if the nodal line of the SC order 
parameter does not cross the FS. Such a case is exactly 
exhibited above in the underdoped regime where the FS 
is lacking around (tt/2, tt/2). Therefore, even if the SC 
order parameter is d-wave which vanishes along the line 
k x = k y , the quasiparticle excitation is still gapped due 
to the absence of the FS across that line. 

As indicated above, it is of general interest to investi- 
gate how the quasiparticle property varies with the in- 
trinsic FS evolution with doping in e-doped cuprates. 
This guides us to conceive the following scenario to 
largely reconcile the experimental discrepancies. Sim- 
ply we assume that the SC order parameter is d-wave at 
all doping levels. This is consistent with the phase sen- 
sitive experimentsii*i£ on optimally doped NCCO and 
PCCO and overdoped NCCO as well as the ARPES 
measurements^ on optimally doped NCCO. In view of 
the FS evolution, we expect that the quasiparticle exci- 
tation will have a finite, i.e., s-wave-like gap in the un- 
derdoped regime and zero gap in the optimally and over- 
doped regime. Then those penetration depth measure- 
ments to support the s-wave gap might be consistently 
understood if nominally optimally or overdoped samples 
in those measurements are in reality underdoped. 

The paper is organized as follows. In the next section 
the model Hamiltonian is introduced and the diagonal- 



ization is performed. In Sec. Ill, the FS evolution from 
underdoping to overdoping, and simultaneously the vari- 
ation of the quasiparticle property and electronic density 
of states (DOS), are exhibited. A remarkable new find- 
ing, i.e., a double-peak feature, in the DOS at optimal 
doping is discovered. In sections IV and V the physi- 
cal quantities heat capacity and superfiuid density are 
calculated, respectively, and their low temperature de- 
pendences are discussed. Conclusive remarks are given 
in Sec. VI. Two appendices A and B are supplemented 
to describe the details of Hamiltonian diagonalization 
and derivation of the electronic spectral function, respec- 
tively. 



II. MODEL HAMILTONIAN AND 
DIAGONALIZATION 

Our idea itself is independent of the concrete model 
calculations as long as the argued FS structures at dif- 
ferent doping levels are recognized. The microscopic 
Hubbard and t-J type models, as extensively applied 
to e-doped cuprates, are difficulliSLS to reproduce the 
global phase diagram. Alternatively, we take as our start- 
ing point the experimental phase diagram and for con- 
venience adopt the phenomenological t-t'-t"-U-V model 
with repulsive on-site U and attractive intersite V to 
simulate the experimental phases. For the F-term 
—V^2(ij) n i n j> only the SC pairing interaction is re- 
tained. The model Hamiltonian reads as follows after 
mean-field decoupling (without irrelevant constant) 

H = -* E (4<*r + H.c.) -t'J2 (clc ja + H.c.) 

(ij)a (ij)'" 

~t" E (4^+H.c.)-]>>m(-l)V + M hU ff 

(ij)"a icr 

-(V/2) ^[A*.(c iT c^ - c a c jT ) + H.c] 

(ij) 

+iV[/m 2 + (V/2)£|A l ,| 2 , (1) 

(ij) 

where (), ()', ()" represent the nearest neighbor (n.n.), 
second n.n., and third n.n. sites, respectively, m = 
(— l)*(ct.Cij — c].Cij_)/2 is the AF order parameter, Ay- = 
(ci-\Cj[ — CiiCji) = 2(cijCjj) is the pairing order parame- 
ter, N is the total number of lattice sites, and the rest of 
the notation is standard. As mentioned above, we assume 
throughout the doping range that the pairing symmetry 
is standard d wave, i.e., A w - = A (—A) for bond (ij) 
along x (y) direction. We set t as the energy unit and 
&b = h = a = 1 (a: lattice constant). Typical values 
t' = —0.3 and t" = 0.2 are adopted. 

To be transparent, the original lattice is divided into 
two sublattices D and E, with net spin orientation ] and 
L respectively. Correspondingly new fermionic operators 
d and e for the two sublattices are introduced, i.e., c$ CT = 
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dia (eia) for i £ D (E). Then the Hamiltonian Q can be and the matrix Hk is 
written in momentum space as 

H = Y^${H k $ k + 2(e' fe - fi)] + N(Um 2 + VA 2 ), (2) 
k 

where the four-component operator ^j. is defined as 

H = (4t e -H 4t d - k i) 
i 



H h 



H — Um 


Afc 







A fc 


-(4- M- 


[/to) 


-Efc 







e' fc — /i + J7m 


Afc 





-£fe 


Afc 


-(4 -M- 



(3) 



with 



4 
Afc 



-2i(cos fc x + cos k y ), 



(4) 



-4i' cos k x cos fc a — 2i"(cos 2k x + cos 2fcj,), (5) 



UA(cos k x — cos k y ). 



(6) 



Throughout the paper the summation with or with- 
out a prime on it means that k is within the magnetic 
Brillouin zone (BZ) — 7r < k x ±k y < tt or the original BZ, 
and the quantities with a hat represent matrices. Under 

a unitary transformation ^ = Uk^k with the new oper- 
ator 

the Hamiltonian (J2J can be diagonalized into 



ff = ^[I'l^fcl'fc + 2(4 - M )] + N{Um 2 + VA 2 ), (7) 



where Hk is a diagonal matrix 
(E7 



Hk — 



-El 



E. 



V 



El J 



with the quasiparticle energy bands 



Ei = J(^-^) 2 + A 2 , 



£ = 4± 



(8) 



(9) 
(10) 



Here the spectra Er show formally the SC pairing in the 
two AF bands which are originated from the begin- 
ning single band in the presence of the AF order. The 
unitary matrix Uk is 



U k = 



COS 0fc 

- sin 4> k 



cos 




sin 


K 


— sin 


% 


cos 




cos 


K 


sinf 




— sin 


K 


cosf 





sm (p k 
cos 4>k 



cos Of sin Of 

-sin 07" cos^ 

cos 0jT sin 6>^ 

— sin (9^ cos 9^ 



(11) 



with </)fc and 9^ satisfying The details of diagonalization are given in Appendix A. 

For late use, we define the 4x4 matrix of the Matsub- 
ara Green's functions 



G(k,r) = -(TAk(T)*l(0)). (16) 



cos20fc = Um/^el + U 2 m 2 , (12) 

sin20fc = Ek/^e 2 + U 2 m 2 , (13) 

cos2# ± — — fi)/E^ (14) Its Fourier transform is easily obtained 

sin2#± = -A fc /£±. (15) G{k,iu n )=U k (iu n -H k )- 1 ul (17) 
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The free energy is given by F = — Tln(Tre~ H / T ) + 
iV(l + x)/i with x: electron doping concentration. This 
leads to 

F = -2T ln[2cosh(^/2T)] 

+N(fix + Um 2 + VA 2 ). (18) 

For given doping x and temperature T, the order param- 
eters m and A are determined by minimizing the free 
energy, where the chemical potential /i is fixed to ensure 
the correct electron number. 



III. FERMI SURFACE EVOLUTION, 
QUASIPARTICLE PROPERTY AND 
ELECTRONIC DENSITY OF STATES 
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FIG. 2: (Color online) (a)-(c) Quasiparticle spectra E k at 
T — along the line T(0, 0)-M(tt, vr) and (d)-(f) Fermi sur- 
faces for the corresponding AF bands g£ — fi plotted in the 
first quadrant [X=(7r, 0)]. In each panel the solid line is con- 
tributed by band "— " and the dashed one by "+". (g)-(i) 
Density plots of the corresponding spectral intensity at T = 0, 
obtained by integrating A(k,uj) [Eq. 1191 1 over an energy in- 
terval [—0.04, 0]t. The broadening factor is O.lt. 

In this work we are focused on the SC (or SC+AF) 
phase and choose three typical doping levels for study: 
x = 0.13, 0.15 and 0.17 which correspond to underdop- 
ing, optimal and overdoping, respectively. In order to re- 
produce the experimental orders at these doping levels, 
the model parameters U and V are slightly changed with 
x: (i) U = 3.9, V = 1 at x = 0.13 (ii) U = 3.8, V = 1.1 



at x = 0.15 (iii) U = 3.7, V = 1.1 at x = 0.17. The 
underlying strategy is to accept part experimental facts 
from the beginning and based on them to work out other 
theoretical results for experimental comparison. The cal- 
culated order parameters m and A at T = are given 
in Fig. ^b), which show the experimental features in- 
cluding (i) the coexistence of both AF and SC orders at 
x = 0.13 and 0.15 (ii) the decreasing AF order with in- 
creasing doping and the vanishing AF order at x = 0.17 
(iii) the strongest SC order (i.e., biggest A) at x = 0.15. 
Accordingly, the quasiparticle spectra EJr at T — are 
plotted along the diagonal line k x = k y (where A& = 0) in 
Figs. Ufa)- (c). It is realized that a minimal finite gap ex- 
ists for x — 0.13 but nodal excitations are present around 
(it/2, 7t/2) for x > 0.15. To thoroughly understand these 
results, we have plotted the Fermi surfaces for the cor- 
responding AF bands £jr — fi in Figs. Htd)-(f). Also, for 
direct comparison with the ARPES data, we have calcu- 
lated the electronic spectral function A(k,u>) (see below) 
and shown the intensity maps in Figs.|2fg)-(i). Explicitly, 
three different cases are distinguished: 

(i) . At underdoping x — 0.13 the FS is lacking around 

(tt/2, 7r/2) [Fig.EJd)]. That is why the finite energy 
(s- wave-like) excitation appears despite the d-wave 
SC pairing. At the same time, it should be empha- 
sized that, due to the proximity of the lower AF 
band to the Fermi level, some spectral intensity is 
already contributed around (tt/2, tt/2) [Fig. |2fg)], 
in consistence with the ARPES dataSS in Ref. |23 

(ii) . At optimal doping x = 0.15 the FS pocket 

around (tt/2, tt/2) truly forms [Fig. Hfe)]. Corre- 
spondingly, the spectral intensity becomes enough 
strong around (tt/2, tt/2) [Fig. Hfh)], as revealed 
by ARPES in Ref. Also note that the FS 

pocket around (tt/2, tt/2) is separated from the oth- 
ers around (7r,0) and (0, it). 

(iii) . At overdoping x = 0.17 a large FS emerges around 

(tt,tt) [Fig. HIT)]. With vanishing AF order, the 
two AF bands will essentially merge into a single 
band and the several FS pieces contributed by the 
two bands will recover to a large continuous curve. 
At present the ARPES data on overdoped samples 
are scarce. According to the evolutionary trend, 
the ARPES intensity map is expected to look like 
Fig. Hi). 

In both cases (ii) and (iii), gapless quasiparticles are in- 
duced from the cross of the FS and the nodal line of the 
SC order parameter (or SC gap) A&. 

Above, the spectral function is calculated by A(k, ui) — 
(-l/TT)lmG(k,iuj n )\ iuJn _ UJ+i0 +. G(k,iuj n ) is the Fourier 
transform of the Matsubara Green's functions G a (k, r) = 
— (T r Cfc (T (T)cJ. cr (0)), which turns out to be spin a- 
irrelevant. For details, one can see Appendix B. A(k,uj) 
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is derived as follows, for all k within the original BZ, 
A(k,oj) = - ^ (l + ^sin20 fe )[(l + cos20£)5(w-££) 

v=+,- 

+(1-cob20(u+£?)]. (19) 



(a)*=0.13 (b)x=0.15 (c)jc=0.17 




FIG. 3: (Color online) Electronic density of states at low 
energy (in unit of t ~ 300 meV) and T — 0. In each panel, 
the red (dark gray) , blue (black) and green (gray) curves show 
the results under the broadening factors V/t = 10~ 4 , 10 -3 and 
5 x 10~ 3 , respectively. 

Furthermore, the electronic DOS per spin is given by 
A^w) = Sfc J 4(^> a -')) which is simplified into 

N(u) = ~ J2 [(l + cos2^)^-^) 

k,u=+,— 

+(1- cos 20£)% + (20) 

The electronic DOS at low energy is plotted in Fig. 
where the results obtained with different broadening fac- 
tors r are given. Let us focus on the curves with 
Y/t = 1CP 4 [red (dark gray) lines]. As expected, the 
DOS takes on a "U" shape at x = 0.13 indicating a fi- 
nite energy gap, whereas it has a "V" shape around zero 
energy for x > 0.15 showing vanishing gap. 

At this stage, we would like to point out another re- 
markable finding in the DOS, i.e., the double-peak struc- 
ture in the case of x = 0.15. (Only the half side u> > 
is talked about.) Careful investigation shows that the 
small peak at lower energy is the coherence peak for 
the "— " band, i.e., the peak position corresponds to the 
SC gap maximum opened on the small FS piece around 
(7r/2,7r/2) [solid line in Fig. [He)], and the large one at 
higher energy is the coherence peak for the "+" band, i.e., 
the peak position corresponds to that on the FS pieces 
shown by the dashed lines in Fig. [2Jc) . Thus this double- 
peak structure is a characteristic feature for the e-doped 
cuprates at optimal doping which have a segmented FS 
and for the two-band modeling of them. In contrast, at 
x = 0.17 the two bands recover to the single band and the 
FS turns into a large continuous curve. Then only a single 
coherence peak shows up [Fig. Etc)], a t the energy cor- 
responding to the SC gap maximum on this whole large 
FS. So far, the double-peak structure has not yet been 



observed in tunneling experiments on optimally doped 
NCCO. 8 ' 9 One possible reason is the effect of the broad- 
ening factor r. As seen from Fig. [21 with increasing Y 
the sharp peaks will be smoothcned and in particular, the 
small peak at the lower energy in the case of x = 0.15 will 
be smeared out. Actually, for Y = 5 x 10~ 3 i ~ 1.5 meV 
which is in the order of the corresponding experimental 
value. 9 the DOS shape at x = 0.15 is more suggestive of a 
s-wave-like gap. Apparently, much smaller Y is required 
in order to observe the double-peak structure and further 
experiments are needed to clarify the issue. 

After extracting the properties of the quasiparticle ex- 
citation and electronic DOS at different doping levels, 
we continue to study their consequences to the physical 
quantities in the following. 

IV. HEAT CAPACITY 

First we calculate the heat capacity, which is given by 
C = T(dS/dT) with the entropy S = -dF/dT. Explic- 
itly S is expressed as 

5 = 2 J2 W cosh(££/2T)] - {EZ/2T) tanh(££/2T)} . 

k, !>=+,— 

(21) 

For general T, the derivative dS/dT has no explicit form 
due to the T-dependence of the order parameters. Only 
at low T <^T C (the SC transition temperature), the or- 
der parameters are nearly T-independent, leading to the 
simple formula for C, i.e., 

C=(1/2T 2 ) J2 (^) 2 /cosh 2 (^/2T). (22) 

k,u=+,- 

The heat capacity per site is shown in Fig. 0] The 
obtained transition temperatures at three doping levels 
are around 30 K, close to the experimental data. That 
the C values at x = 0.15 are bigger than those at x = 0.17 
in the region about T/t < 0.002 is consistent with the fact 
seen in Fig. [31 that the DOS at x = 0.15 is greater than 
that at x = 0.17 in the very low energy region [below the 
small peak in Fig. 13Tb)]. Here we are mainly interested in 
the qualitative behaviors at low T <C T c . It is found that 
the heat capacity takes on an exponential T-dependence 
for x = 0.13 due to the s-wave-like excitations, but a 
power law behavior for x = 0.15 and 0.17 due to the 
nodal excitations. To more clearly see this point, we 
have replotted the low-T data for x = 0.13 in the lower 
inset by using — 1/ \n{C /Nks) vs T/T c and those for x = 
0.15 and 0.17 in the upper inset with logarithmic scales 
for both axes. Now we have nearly straight lines in all 
cases. Moreover, it is easy to check that the slopes for 
the two lines in the upper inset are both 2, which means 
T 2 dependence of C for x = 0.15 and 0.17. In addition, 
we find that in the pure AF phase without SC order, e.g., 
at doping x = 0.1, the heat capacity takes on a linear In- 
dependence (not shown), indicating an AF metal. So it is 
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Pk = UlhUk, (28) 



gk = -d 2 £k/dk1 and g' k = —d 2 e' k /dk x . Then the ther- 
mal average is obtained 



(jfx) = EEA)»/(i+^) 



(29) 
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FIG. 4: Specific heat as a function of temperature at three 
doping levels. It is seen T c /t = 0.008, 0.0113 and 0.0075 
at x = 0.13, 0.15 and 0.17, respectively. The data at low 
temperatures (T/T c <C 1) are replotted in the lower inset for 
x — 0.13 and upper inset for x — 0.15 and 0.17. Note that 
the different axis rescales are used in the two insets. 



of great interest to experimentally extract the electronic 
part of the heat capacity and investigate the variation of 
its T-dependence with doping. 



where (Pk)ss' represents the matrix element of Pk and 
for convenient statement the quasiparticle energy bands 
are relabeled by E% (s = 1, 2, 3, 4) with E\ = E k , E\ = 
-E% , El = Et and Ef = -E+. 



down the current operator , i.e 



To calculate the second term of Eq. I|23|l. we first write 
; operator J*, i.e., 



(30) 



with 



V. SUPERFLUID DENSITY 

Another important quantity for calculation is the su- 
perfluid density p which relates to the penetration depth 
A by p cx A~ 2 . The superfluid density divided by the 
carrier effective mass m* is given bj^S 

p/m* = (-K x )/N - A(q x =0,q y ^ 0,iO„ = 0), (23) 

where the operator K x is derived as follows (x, y: unit 
vectors in x, y directions, respectively) 

K x = -t£(4 fa ,ci a +H.c.)-t' £ + Hx -) 

-^"E^W^+H.c.) (24) 

la 

and the current-current correlation function A is (/3 = 
1/T and fl n = 2nir/[3: bosonic Matsubara frequencies) 

A(q,iQ n ) = ^ j\r^{T T J-{T)r_ q {Q)) (25) 

with J% = e~ lqRl Jf. Here Jf is the x-component of 
the paramagnetic current operator at site I 

J f = ^£(4+^* - H - c -) + **' £ (4+5v c ^ - H - c -) 



Vk, q = 





fc 


^ f'k,q 








- f 

J —k+q,q 


fk,q 





V 


~ f—k+q,q 



fk,q 
-f-k+q,q 







o -/: 



/fc,. 



j t [ e -»(fc*-g*) _ e 



(31) 



and 



Above, 

f k q = 2it'[e~ l ^~ q ^ cos(k y - q y ) - e lk * cos k Vj 

2it"[e~ 2l ( k *~ q ^ — e 2ifex ], which reduce to the simple 
form at q = 0: / fej0 = de k /dk x and /£ = de' k /dk x . 
According to Eq. I|25[) . the correlation function A is 
derived as, after expansion of G{k, t) into G(k,iu} n ), 



A(q,m n ) = -^££ Tr [G(fc ~ 1^ iuJ n)Vk,q 
k iu> n 

G(k,iQ n + iu} n )Vk-q,-q\- (32) 



Substituting G(k, iu) n ) with Eq. (|17(l and completing the 
summation over iu> n , we obtain 



1 - 
A(q,ifl n ) = t^E E (^,9)^(^-9 

fe s,s ; — 1 

tanh(£|_ ./2T)-tanh(^/2T) 



iV n + - E( 



,(33) 



-2^"£(4 +25(7 c i<r -H.c.). 



(26) 



The first term of Eq. 12: il) can be easily calculated if 
K x is expanded in momentum space, i.e., 



where (Wk.q)ss' represents the matrix element of Wfc, q = 
Ul_ q Vk. q Uk- Finally we set iCl n = and then take the 
limit q y — > to have 



K * = E H p ^k = E hPk*k, 
k k 



(27) 
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A(g x = 0, % -> 0, 4 O n = 0) = — 22 \z2 K W k,o)ss>] 2 fc ' _ ps > 



[(Wkfi)s 



2N kti{^ s E k~ E k 2Tcosh^/2T) 



m 




T/T 

c 



FIG. 5: Superfluid density as a function of temperature at 
three doping levels. The low-T regime is enlarged in the inset. 
The circles in the main panel are the experimental data for 
PCCO at x — 0.131 extracted from Ref. (see text). 

The superfluid density is given in Fig. |5j where the 
ratio p(T)/p(0) vs T/T c is plotted. Still we are concen- 
trated on the low-T regime, which is enlarged in the inset. 
The notable observation is that the curve for x = 0.13 
is very flat, which may be characterized by the formula 
1 - p(T)/p(0) oc exp(-6/T) with a constant b. When 
the quasiparticle excitations become gapless at x = 0.15 
and 0.17, the results are much different. Now p(T) de- 
creases faster, with a linear relation 1 — p(T)/p(0) oc T, 
as apparently seen from x — 0.17 in the inset of Fig. 
At x = 0.15, p(T) / p(0) in the main panel shows an in- 
flexion point as a result of crossover from x — 0.13 to 
x = 0.17. Similar to the case for x — 0.17, the linear 
relation holds at T/T c < 0.1 due to the nodal excitations 
on the small FS pocket around (ir/2,ir/2) [solid line in 
Fig. life)]. When T increases the low energy quasiparti- 
cles around this small FS pocket are completely excited. 
Then the change of the superfluid density is controlled by 
the finite energy excitations around the FS pockets cir- 
cling (ir, 0) and (0, 7r), similar to the case for x = 0.13. So 
the curve at x — 0.15 away from low-T tends to behave 
like the one at x = 0.13. The linear low-T-dependence 
at x — 0.15 and 0.17 may become T 2 in real materials if 
impurity scattering is considered^ 

The low-T-dependence of the superfluid density (or 
penetration depth) is important to hint the SC pairing 
symmetry and has been substantially measured for e- 
doped cuprates. 6 ' 7 - 15 ' 16 ' 20 However, the experimental re- 
sults are controversial and some of them might be mi- 



sexplaincd. For example, the s-wave SC gap was in- 
ferred, from the flat feature of the experimental data, 
by Kim el al£ for PCCO over a wide doping range 
0.124 < x < 0.144. In Fig. our theoretical result at 
x = 0.13 is compared with their data at x = 0.131, and 
good agreement is obtained. But actually in our theory 
the SC gap Aj, is d-wave. It should be reminded that the 
doping level x = 0.131 in the work by Kim et al. cor- 
responds to the optimal doping, although this optimal 
value is 0.15 for PCCO in most experiments4^*ii*i^ Our 
result suggests that, even with a d-wave SC pairing, it is 
not difficult to understand the s-wave- like excitation from 
some penetration depth measurements at various doping 
levelsSiSS if nominally optimally or overdoped samples 
in those measurements are in reality underdoped, or more 
generally, if the FS is absent around (tt/2,it/2) in those 
samples perhaps due to the still strong AF order. On the 
other hand, the power law T-dependence (T or T 2 ) of p at 
optimal and overdoping is consistent with other penetra- 
tion depth measurements^ 5 ^ and many others such as 
phase sensitive Josephson junctionsii*^ and ARPES. 13 
Thus in this way the diverse experimental results could 
be largely reconciled. Finally, we mention that our es- 
sential idea to exhibit the s-wave-like excitation and to 
explain the flat feature of some penetration depth data is 
different from the phenomenological superfluid response 
theory in Ref. where the existence of nodal quasipar- 
ticle excitations is assumed at all doping levels. 



VI. CONCLUSION 

In conclusion, we have explored an intrinsic physical 
mechanism, based on the FS evolution with doping, to 
understand the diverse experimental results on the SC 
pairing symmetry in e-doped cuprates. We argue that 
the FS pocket around (7r/2, 7t/2) has not yet formed un- 
til doping reaches about the optimal value. Therefore, 
in the underdoped regime, even if the SC pairing order 
parameter is d-wave which vanishes along the diagonal 
line, the quasiparticle excitation is still gapped which 
shows a s-wave-like behavior due to the absence of the 
FS across that line. To demonstrate this idea, we have 
elaborated a theory with consideration of coexisting AF 
and SC orders in e-doped cuprates. The variations of the 
quasiparticle property and electronic DOS are exhibited 
as the FS evolves from underdoping to overdoping. The 
physical quantities heat capacity and superfluid density 
are calculated and their low temperature dependences 
are compared at different doping levels. Our study gives 
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new understanding of some measurements which actually 
detected the s-wavc character of the quasiparticle exci- 
tation (rather than the SC order parameter) and thus 
tends to reconcile the apparently contradictory experi- 
mental results. We have also predicted the remarkable 
double-peak feature in the DOS for the e-doped cuprates 
at optimal doping which have a segmented FS. 
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APPENDIX A: HAMILTONIAN 
DIAGONALIZATION 

We consider the following two steps to diagonalize the 
Hamiltonian J3J, which is explicitly 



Second, we adopt the Bogliubov transformations 



H = Y, £ k( d l^ka +H.c.)+J](e , fc -/i~aC/m)4 (T d 

ka ka 
I 



k 

+N(Um 2 + VA 2 ). 



(Al) 



First, we use the rotational transformations to remove 
the hopping term between d and e, i.e., 



' ] H, ( t k ] ), f T. )=U lk ( |j ) ,A2, 



efet J ' V dk l 
with [see Eqs. (I12[l and i|13|) for the choice of 4> k ] 

U lk =( 008 1 Shi t k ). (A3) 
^ -sm^fe cos^fe J y ' 

In terms of the new operators d and e, (| Al|) turns into 

k 

k 

- ^2 A fe (J feT e_ fci + e k -\ d-ki + H.c.) 

k 

+N{Um 2 + VA 2 ), (A4) 
where the pairing term is formally unchanged. 



'2k Jt 



-kl 



with [see Eqs. (1141) and 115|) for the choice of i 



(A5) 



(A6) 



Then l|A4|l is diagonalized into 

H = J2^ E k(^ k1 dk]- - e- kl el kl ) 

k 

+E+(el r e m - d- H $_ k] ) + 2(4 - »)] 
+N(Um 2 + VA 2 ), (A7) 

which is exactly the Hamiltonian Q. If all the transfor- 
mations l|A2(l and (| A5|) are integrated, we will arrive at 

the global transformation ^Sf k = U k ^ k with U k given by 
Eq. (TTJ. 



APPENDIX B: DERIVATION OF THE 
SPECTRAL FUNCTION EQ. (tT51l 

With division of the original lattice into two sublat- 
tices, one has the following correspondence for the oper- 
ators c and d/e in momentum space [Q — (tt,tt)] 



Cka = {dka + eka)/V2, 
Ck+Qa = {d ka - e,ka) / 'V%, 



(Bl) 
(B2) 



where k is restricted into the magnetic BZ. Then the 
Green's functions for the original c-electron can be ex- 
pressed as follows, with help of the matrix elements [de- 
noted by G ss ' (k, t)] of G(k,r), 



G t(^) = ~ £ G ss ,(k,r), 



(B3) 



s,s / = l,3 



G.i(k,r) = - 1 - G s A~k,-r), (B4) 

G T (fc + Q,r) = ~ tss'G ss ,(k,T), (B5) 

s,s'=l,3 

G;(fc + Q,r) = ~ Yl *ss'GsA-k,-T), (B6) 

s,s'=2,4 

where the coefficient e ss ' — 1 (—1) for s' = s (s' ^ s) 
is introduced. All the above equations are also valid for 
the Fourier transformed functions after the replacement 

T — * lLO n . 

By use of Eq. (|17fl . we can obtain the following sums 
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h(k,iu n ) = G ss ,{k,iu n ) = J2 (l + ^sin20 fc ) ( . cua ~ P *„ + . Sm ^ ) , (B7) 

I 2 (MO = J2 G s Ak,iw n ) = V (1 + v sin 2<f> k ) ( . °l v + . C0S ' „ ) , (B8) 

s,s'=2,4 v=+,- \ ™ fc « 1 fc / 

7 3 (MO = J £ ^ tf W= + (B9) 

h(k,iu n ) = J2 e ss ,G ss ,(k,iaj n )= V (l-i/sin2</. fe )f , Sm g * „ + • ^ V ( B1 °) 

, n „ ^ lUJn+EYJ 

s,s'=2,4 1/= + ,— x K « * 

It is easy to check Gj(fc,iw„) = — |i2( — &j —iu n ) — —^l2(k,—iui n ) = \lt{k,iw n ) = Gf{k, icu n ). Similarly one has 
G^(fc + Q, iuj n ) — Gj(fc + Q, io; n ). So the spin index can be dropped out and we obtain 

G(k,iuj n ) = ij^fcjiu^), (Bll) 
G(k + Q,iuj n ) = ^I 3 (k,iu n ) = ^I 1 (k + Q,iu} n ), (B12) 



r 



where the last equation holds if the definitions for 4>ki G(k,iu> n ) has the unified form i|BIl|) . Analytical contin- 
6fr and Efr are extended to the original BZ. There- uation iui„ — > + z0 + for Gffc. ia; w ) will lead to Eg. {T5J| . 
fore, no matter fc is within or beyond the magnetic BZ, 
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